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Abstract 

We present an experimental study of velocity statistics for a partial layer of 
inelastic colliding beads driven by a vertically oscillating boundary. Over a 
wide range of parameters (accelerations 3-8 times the gravitational accelera- 
tion), the probability distribution P(v) deviates measurably from a Gaus- 
sian for the two horizontal velocity components. It can be described by 
P(v) ~ exp(— |v/u c | 15 ), in agreement with a recent theory. The characteristic 
velocity v c is proportional to the peak velocity of the boundary. The granular 
temperature, defined as the mean square particle velocity, varies with parti- 
cle density and exhibits a maximum at intermediate densities. On the other 
hand, for free cooling in the absence of excitation, we find an exponential 
velocity distribution. Finally, we examine the sharing of energy between par- 
ticles of different mass. The more massive particles are found to have greater 
kinetic energy. 

PACS: 83.70.Fn, 05.20.Dd, 45.40.-f, 83.10.Pp 
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We determine the statistical properties of particles in a vibrated granular medium ex- 
perimentally. While many similarities between ordinary gases and excited granular media 
have been found, a fundamental difference is that collisions between particles in granular 
matter are inelastic. As a consequence, the velocity distribution deviates measurably from 
a Gaussian, but can be described by P{v) ~ exp(— |f /fc] 1 ' 5 ) for a large range of parameters 
where external excitation is sufficiently frequent, in agreement with a recent theory. 



I. INTRODUCTION 



Granular systems are often treated statistically, since the large number of degrees of 
freedom and the complexity of interparticle forces limit analyses on the particle level. What 
are the statistical properties of particles in an excited granular medium? In many respects 
the dynamics of granular media strongly resemble the dynamics of ordinary fluids and solids. 
Yet there are fundamental differences between granular solids and fluids and their molecular 
counterparts that make the answer to this question intriguing 0,0. Ordinary temperature 
does not lead to measurable velocity fluctuations due to the large mass of granular particles. 
When granular particles are excited by an external energy source, their thermal energies are 
much smaller than their kinetic energies. During collisions between particles, frictional forces 
and deformation near the points of contact lead to dissipation of kinetic energy. Therefore 
it is not surprising to find that kinetic energy is transformed into thermal energy when 
particles interact in an excited granular medium [Q. 

For strongly excited granular media, assumptions similar to those in kinetic theory are 
often made 0. Among these assumptions are that velocity distributions are Gaussian and 
that the mean energy is shared equally among the various degrees of freedom. However, 
recent numerical and theoretical research indicates that excited inelastic hard spheres can 
exhibit non-Gaussian velocity distributions PHlQf, though the predicted velocity distribu- 
tions differ considerably from each other. Some thermodynamic descriptions of granular 
media make use of the concept of entropy |TT[, or separate the dissipative degrees of free- 
dom from conservative ones . 

Vibrated granular media have played a special role in efforts to understand the dynamics 
of granular materials, in part because vibration is a convenient method of replacing the 
energy lost to friction and inelasticity. A variety of novel phenomena have been discovered 
in the past fifteen years [0: heaping and convection rolls []l3| , |l^| , standing and traveling 

waves 



T^j-fT?!, oscillons and fluidization |I5| - |2"I]| . 

In this paper, we study the simple case of a single layer of particles with vertical excita- 
tion, varying the fractional coverage of this layer and the amplitude of the driving waveform. 
We determine the velocity distributions precisely and compare them to various functional 
forms, an issue that previous experiments have left unresolved for granular particles that are 
free to move in three dimensions. Quantitatively studying particles in three dimensions has 
only recently been achieved in rotated cylinders using NMR techniques; it yielded important 
information about the segregation process p2| |. By observing a three dimensional excited 



granular medium from above, we are able to focus on the shapes of the distributions of 
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the horizontal velocity components. A combination of otherwise identical white and black 
particles allows us to track individual particles even for high fractional coverages at high 
excitations where particles frequently move over each other. We measure the dependence of 
the variances of the distributions on the excitation amplitude, frequency and the fractional 
coverage. Sharing of energy in mixtures of granular materials is also investigated. 

In Sec. II, we discuss the background for our experiments. In Sec. Ill we present the 
experimental setup. In Sec. IV, we show that, in most cases we treated experimentally, 
the velocity distribution P(v) deviates measurably from a Gaussian and can be described 
by P(v) ~ exp(— |f /f c | 15 ). We show precisely how the variance increases with driving 
amplitude and changes with fractional coverage c, where c = 1 is the coverage for a compact 
crystal of beads. The energy sharing between particles of different types is also described 
here. 



II. BACKGROUND 



A number of studies on excited granular media have addressed the issues of clustering 
and energy sharing. Numerical simulations and theoretical derivations have shown that 
the presence of inelasticity in granular flows can lead to the formation of clusters; as a 
consequence, equipartition of energy fails [p3|-p5|. Experiments performed in a horizontal 



two-dimensional layer were consistent with this predicted clustering effect |26 |. Vertical 



one-dimensional experiments and simulations were performed earlier |2U||; a crossover from a 



condensed (clustered) to a fluidized state was found as a function of the driving acceleration, 
the number of beads, and the coefficient of restitution. 

Deviations from equipartition due to clustering are straightforward to understand on 
physical grounds. Inelastic collisions imply a loss of energy each time a collision occurs. 
When particles begin to gather in a certain region of space, the rate of their collisions 
increases. The rate of energy loss for this group of clustered particles is thus greater, and 
the distribution of their velocities becomes narrower than that of the particles in less dense 
regions. 

Clustering can affect the velocity statistics. Assuming a Gaussian velocity distribution 
for a nearly homogeneous granular medium, Puglisi et al. f7||| predict a non-Gaussian ve- 
locity distribution due to clustering: a superposition of Gaussian velocity distributions with 
different widths. For inelastic particles this leads to high velocity tails in the velocity distri- 
bution which decrease more slowly than a Gaussian function but faster than an exponential. 
The velocity distributions obtained in a simulation by Taguchi and Takayasu |27| are power 



laws, resulting from clustering. These distributions have diverging variance; this calls into 
question the notion of a granular temperature, which is generally derived from a variance. 
In a two-dimensional simulation, Peng and Ohta [[J found that the velocity distributions 
deviate from Gaussian behavior under the influence of gravity unless g8h << T, where Sh 
is the height of the region of observation. T is the granular temperature (see e.g. Ref. JT|), 
defined as the variance of the velocity (minus the mean velocity) 
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T =< (v 2 - <v> 2 ) > 



(1) 



The high energy tails of velocity distributions in a homogeneous granular fluid were 
investigated theoretically by Esipov and Poschel || for the unforced case, and by Noije 
and Ernst ]6[] for both the unforced and the heated case, based on the Enskog-Boltzmann 
equation. In the free case, i.e. without energy input into the system, the velocity distribution 
approaches an exponential at high velocities. When energy is added to the granular medium 
randomly and uniformly throughout the system, the high velocity approximation becomes 
P(v) ~ exp(— \v I 1 ' 5 ). In order to compare our experimental results to this prediction we 
have to assume that energy input into the measured horizontal motion occurs randomly. 

While measurements of velocity statistics in an excited granular medium have been 
carried out |28| , few measurements precise enough to distinguish between different functional 
forms of the velocity distribution exist to our knowledge. One exception is a recent study 
of clustering and ordering near a peak acceleration of a = 1 g by Olafsen and Urbach |29| . 



who found significant deviations from Gaussian distributions, at low and especially at high 
velocities, where the distributions become exponential. 

Deviations from equipartition can occur for reasons other than clustering. For exam- 
ple, Knight and Woodcock PD[ studied a vibrationally excited granular system theoretically 



(without gravity) and concluded that equipartition need not be observed at high amplitude 
of excitation due to the anisotropy of the energy source. A two dimensional system (a verti- 
cal Hele Shaw cell) was studied experimentally by Warr, Huntley, and Jacques |28]. Velocity 
distributions, though roughly Gaussian, were checked and found to exhibit anisotropy be- 
tween the vertical and horizontal motion: the horizontal velocity distribution was narrower 
than the vertical one. Grossman, Zhou, and Ben-Nairn |3l[] considered a two-dimensional 



granular gas (without gravity) from a theoretical quasi-continuum point of view and found 
the density to be nonuniform and the velocity distributions to be asymmetric for "thermal" 
energy input from one side. McNamara and Luding |53] considered the sharing of energy be- 
tween rotational and translational motion of the particles, and found a significant violation 
of equipartition. 

The scaling of the granular kinetic energy with vibration amplitude is also of interest in 
connection with the experiments to be discussed in the present paper. It has been considered 
experimentally by Warr et al. [28|, numerically by Luding, Hermann, and Blumen and 



theoretically by Kumaran in the nearly elastic limit of weak dissipation |3|] and also by 
Huntley in a simple model. The results of these different studies do not seem to be 
mutually consistent with each other, perhaps because different regimes were explored; the 
situation is unclear. 



III. EXPERIMENTAL SETUP AND METHODS 



The experiments are conducted in a circular container of diameter 32 cm made of delrin. 
It is driven vertically with sinusoidal acceleration at a single frequency using a VTS500 
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vibrator from Vibration Test Systems Inc. A computer controlled feedback loop keeps the 
vibration amplitude constant and reproducible. The frequency / used in most experiments 
is 100 Hz and the peak acceleration a of the plate is in the range 3 — 8 g. The peak plate 
velocity v p = a/2-nf lies between 4.7 cm/s and 14.1 cm/s. The particles are glass beads 4 
mm in diameter (from Jaygo Inc.) with fractional coverage c. A glass cover at 2 cm height 
prevents the beads from escaping from the container. Collisions with the cover are rare for 
most experimental parameters but lead to measureable changes in the velocity distributions 
at the largest accelerations, if the coverage is low. Some charging of the glass beads is 
noticeable when the beads are at rest, but electrostatic forces are negligible at the range of 
accelerations investigated here. 

Our objective is to measure the horizontal velocity distributions in an excited three di- 
mensional granular medium for a large range of particle densities. However, particle tracking 
becomes increasingly difficult as particle tracks approach each other and cross frequently at 
large c, thus hindering reliable identification of the horizontal positions of individual parti- 
cles. However, if only a modest number of particles are reflective, the frequency of collisions 
between them is small and nearly independent of coverage or excitation. By tracking these 
test particles we can also directly compare the velocity distributions for different coverage. 

Tracking a subset of the particles is accomplished by using some white glass beads among 
black glass beads. Except for the color, the black and white beads have identical physical 
properties. Stainless steel beads of three different diameters replaced the white glass beads 
in some experiments. The material properties of all particles we used are listed in table |I[ 

Images of an area 16.74 x 15.70 cm at the center of the container are taken at a res- 
olution of 512 x 480 pixels using a fast camera (SR-500, Kodak Inc.) operated at 250 or 
500 frames/s. At a vibration frequency of 100 Hz this ensures that images are taken at 
5 different phases relative to the phase of the plate vibration yielding the average energy 
throughout the cycle. The images are analyzed using IDL (Research Systems Inc.) software. 
First each image is enhanced using a bandpass filter and thresholding to eliminate noise. The 
positions of all bright particles are found from the enhanced image by calculation of their 
centroid; this defines the particle positions reproducibly to within less than 0.1 pixels. Small 
effects due to the finite pixel size are noticeable sometimes as slightly increased probabilities 
of particle displacements that are multiples of the pixel width (0.0327 cm in the physical 
system). This displacement corresponds to v — 8.18 cm/s in most experiments. The high 
frame rate ensures that even the fastest beads move less than one particle diameter between 
images. This allows accurate tracking of all bright beads for all 546 sequential frames (the 
maximum available with our camera) with a typical precision of ±2% for the velocity mea- 
surements. We project each step onto two perpendicular directions, and study the statistics 
of each velocity component. One concern was that as a black bead moves over a bright 
bead, the centroid position could move away from the center of the bead and possibly alter 
the measured velocity distribution. For a test we thus eliminated particles whose integrated 
greyscale intensity changed rapidly. We found that eliminating these points does not mea- 
surably alter the velocity distribution. However, if a collision occurs between frames, our 
measurements indicate the average of the velocity prior to and after a collision. The mea- 
sured distribution of velocities will therefore probably be slightly closer to a Gaussian than 
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the real distribution of velocities. To improve the data in this regard, it would be necessary 
to sample faster while retaining the same relative accuracy in the velocity measurement. It 
would therefore be necessary to use a faster camera, to zoom in closer to the sample and 
to extend the measurement over a significantly larger number of frames to obtain the same 
statistics. 



IV. EXPERIMENTAL RESULTS 

A. Granular Temperature 

Extracted particle tracks are shown in Fig. |l] for c = 0.28 (half of the beads are white) 
and a = 5 g. For clarity, only tracks longer than 200 images are shown, which eliminates 
some tracks close to the edge. A total of 546 frames is acquired at 250 frames/s (i.e. for 
2.18 seconds) with approximately 200 tracked particles in each frame. 

The particle velocities are determined from the particle displacement between consecutive 
frames. This does not always represent the true velocity of the particle though. If a collision 
occurs between frames, the apparent velocity will be lower than the true velocity. We can 
define an apparent temperature based on displacements along either horizontal coordinate 
denoted here as x; it depends on the time between frames At: 

T(At) =< ( Xj (t k + At) - xjih)) 2 > j>k /At 2 , (2) 

where the average is taken over all particles (j) and frames (k). The particle tracks 
obtained at 250 frames/s during an interval of 2.18 s allow us to determine T(At) for 
1/250 s < At << 2.18 s. It is often useful to express the velocities in units of the peak plate 
velocity v p , which yields a dimensionless temperature T(At) = T(At)/v 2 . Fig. |2| shows the 
dimensionless temperature T(At) vs. frame rate (1/At) at c = 0.42 and / = 100 Hz for 
different accelerations. For small frame rates, i.e. large At, T(At) increases approximately 
linearly with the frame rate T(At) ~ 1/At. This indicates that the particle motion may be 
described by an ordinary diffusion law when many collisions occur in the sampling interval 
At. Assuming such a diffusion process, one expects that 

T(At) « (for At » r c ) . (3) 

The dashed and the solid lines in Fig. |2| are linear in 1/At and give upper and lower limits 
to Tt c . We estimate that T « 1.0 based on the high frame rate limit. We can now estimate 
the collision time to be 0.02 s < r c < 0.04 s for / = 100 Hz and c = 0.42, which corresponds 
to one collision every 2 — 4 oscillations of the vibrator. On the other hand, for large frame 
rates (small At), T(At) approaches a constant. This occurs when At is much smaller than 
r c ; in this limit T(At) rs T, so the granular temperature T (for displacements along one 
axis) can be defined as 
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T = lim T(At) =< v 2 > 

Ai-+0 



(4) 



The approach of the measured granular temperature to a constant cannot be fitted by 
a simple exponential or power law. For an excited granular material, this shape could be 
influenced by correlations in velocity between neighboring particles and correlations between 
the local density and particle velocity. Starting from the fact that the sampling time is 
roughly 10% of the mean collision time at an intermediate coverage, we estimate that the 
true granular temperature T might be up to 10% higher than the measured value. 

All values of T and T presented in this paper are the measured temperatures obtained 
at the highest frame rate (usually 250 frames/s), where the limit At — > is justified. This 
limit is approached in a very similar way for different accelerations (see Fig. |2|). This allows 
us to determine the acceleration dependence of T. Note that the collision time can be 
roughly independent of acceleration since the height of the average bounce increases with 
increasing peak plate acceleration. For the fixed number of particles in our system this 
leads to an increase in the mean free path with increasing acceleration. When the number 
of particles is changed, the collision time also changes. As shown in Fig. |3], the limit of 
Eqn. (4) is approached fastest at the lowest energy of c = 0.14 and slower at c = 0.42 and 
at c = 0.98. However, we are close enough to the limit at all coverages (since the mean time 
between collisions is always significantly longer than the time between frames) to observe 
qualitatively how T changes with c. 

We can therefore determine the dependence of the temperature on plate acceleration a 
and coverage c, shown in Fig. |j. The temperature increases with a for all coverages. As a 
function of c, T increases at low coverage, exhibits a maximum around c = 0.30 and decreases 
with c at high coverage. This trend reflects changes in the true granular temperature as 
shown in Fig. 0. No measurable change in the temperature dependence occurs around 
c = 1. A fractional coverage above unity is meant to indicate that more than one close 
packed layer of beads is used. The highest granular temperature indicates that the average 
potential energy of the particles corresponds to a mean height above the plate comparable 
to one particle diameter d. All experimental results are therefore limited to the regime of 
particle energies smaller than or comparable to the only characteristic energy of a granular 
material, the potential energy of raising one particle by one diameter mgd. When scaled by 
the peak plate velocity v p as in Fig. |5j, the granular temperature becomes independent of 
acceleration to within ±10% for most data points. Remarkably, the dependence on coverage 
follows approximately the same behavior at all accelerations. 

The scaling of velocities by v p implies that T ~ l// 2 . Fig. |B| shows, on a log-log plot, 
that T does indeed decrease approximately ~ 1/ f 2 for the two accelerations and the three 
coverages shown. The plot covers two orders in magnitude of the granular temperature, 
ranging from conditions where beads rarely hit the container lid to conditions where frequent 
collisions with the lid occur. We conclude that the scaling of the bead velocity by the plate 
velocity is very robust and is not significantly influenced by additional contacts with the 
container lid. 
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B. Velocity Distributions 



The velocity distribution along one axis, obtained from particle tracks of white beads, is 
shown in Fig. ^. Fits to a Gaussian distribution F g (v) = A g exp(— \v/y2v c \ 2 ) are shown as 
dashed lines, and fits to the prediction of Ref. |J 

F 2 (v) = A 2 exp(-\v/l.lQ4v c \ 1 - 5 ) (5) 

are shown as solid lines. The characteristic velocitiy v c is defined as the square root of the 
variance, so that for both F g and F 2 . 

Vc = V< v 2 > = VT. (6) 

The data points for the fits are weighted equally on a linear scale in Fig. |7](a,b) and equally 
on a logarithmic scale in Fig. 0(c,d). The characteristic velocity v c in (a,b) is v c = 6.19 cm/s 
for F g and v c = 6.99 cm/s for F 2 . The Gaussian fit underestimates the probability of both 
low and high (v > 3v c ) velocities, while the fit to F 2 describes the probability distribution 
quite well over three orders of magnitude in probability. The increased weight of the high 
velocity experimental data in (c,d) leads to v c = 7.66 cm/s for F g and v c = 6.81 cm/s for 
F 2 . The Gaussian fit again underestimates high and low velocity probabilities, while the fit 
to F 2 proves to be insensitive to the weighting of data points, indicating a robust fit. The fit 
to F 2 is also insensitive to the choice of the fitted range of velocities, while the characteristic 
velocity decreases with decreasing fitting range for F g . We conclude that F 2 provides a 
better fit than F g . 

The velocity distribution for a large range of accelerations 3 g < a < 8 g can be described 
accurately by F 2 as shown in Fig. ||(a). The data and the best fitting lines are shifted 
vertically as needed for clarity; this amounts to multiplication by a constant on a log-linear 
plot. The probabilities are plotted against |w| 1,5 , where v = v/v p . On this log-linear scale, 
F 2 is a straight line, in good agreement with the experimental data for a large range of 
accelerations. The dependence of the velocity distribution on coverage is more complex. 
The coverage dependence is shown in Fig. |8](b). At relatively high coverage, above c = 0.28, 
the experimental data can be described well by F 2 (solid lines). However, at lower coverage 
(for approximately the same range of c where T increases with c) the probability of high 
velocities is underestimated by F 2 (dashed line). 

While the granular temperature is proportional to 1/ f 2 to a good approximation, mea- 
sureable differences are apparent in the distribution of non-dimensional velocities v, as a 
function of frequency. Fig. ^| compares the distributions for 40 Hz and 140 Hz. The dis- 
tribution falls more slowly with velocity as / is increased. The cause of this behavior is 
probably the smaller (unsealed) velocities at higher frequency, which decreases the collision 
rate (i.e. the rate of energy loss through inelastic collisions), and the higher frequency of 
vibration (i.e. roughly the rate of energy input). For the highest unsealed particles velocities 
and lowest excitation frequency (a = 5 g and / = 40 Hz) deviations from fits to F 2 become 
observable. 
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At low accelerations a = 2 g we observe that the velocity distribution has exponential 
tails and an approximately Gaussian central component. The crossover from a Gaussian 



distribution to an exponential distribution is shown in Fig. [TO]. These results are similar to 
the velocity distributions around a = 1 g obtained by Olafsen and Urbach ||29|| , and possibly 
related to clustering effects observed at low accelerations. In addition, the system is nearly 
two dimensional since most beads bounce lower than one particle diameter.. 

It is possible to investigate the case of a freely cooling granular medium to some extent 
by looking at a time averaged velocity distribution. We shut off the vibrator abruptly, and 
simultaneously trigger the camera at 250 frames/s. We then extract particle tracks from 
the image sequence in the same way as we did for image sequences of continuously excited 
granular media. In order to analyse the functional form of the velocity distribution, we need 



to accumulate velocities over 150 frames, i.e. 0.6 s. Fig. |TT| shows such a velocity distribution 
for c = 0.84, where the vibrator was operated at / = 100 Hz and a = 5 g prior to shut 
off. The velocity distribution is exponential, in agreement with calculations for the free 
cooling case ]5|||. However, the measured velocity distribution represents an average over 
almost the entire free cooling process, since the instantaneous velocities indicate that the 
granular temperature decreases by more than one order of magnitude during the averaging 
time. We have also made measurements over 0.1 s; although the statistics are not as good, 
the distributions still appear to be exponential. 



C. Equipartition 



In a final set of experiments, the white glass beads are replaced by steel beads of different 
size (grade 100 stainless steel 316 for the two smaller sizes, grade 100 stainless steel 302 for 
the largest beads). At a total coverage c = 0.42, 14% of the glass beads were replaced 
by steel beads (for the smallest bead size only 5% of the glass beads were replaced). In all 
experiments most collisions of steel beads therefore occur with glass beads. The distributions 
of v for the tracked particles are shown in Fig. 12. Larger beads have smaller non-dimensional 
characteristic velocities v c and thus a smaller granular temperature T than smaller beads. 
The velocity distributions for the two smaller steel bead sizes are described well by F 2 , while 
that for the largest beads is better described by a Gaussian (dashed line). This could indicate 
that the large beads effectively prevent clustering since they act as a source of momentum 
for the other beads. While the granular temperature of the largest steel beads is lower than 
the temperature of the glass beads, their energy mT is larger, while the energy of the smaller 
steel bead sizes is smaller, as shown in Fig. pj . 



V. SUMMARY AND CONCLUSION 



We have reported experimental studies of velocity statistics for a (fractional) layer of 
glass beads subjected to vertical vibration. The horizontal motion of a small subset of 
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beads was measured using a high speed camera at a frame rate sufficiently high to mea- 
sure instantaneous velocities accurately. The measurements were acquired over an interval 
substantially longer than the time between interparticle collisions. These capabilities al- 
lowed us to determine particle statistics of the horizontal motion in detail. We analyzed 
granular temperatures and velocity distributions for a large range of excitation frequencies, 
amplitudes, and coverages. 

The variance of the particle velocity distribution (or the granular temperature T of 
horizontal motion) varies approximately in proportion to the plate velocity (Fig. |5]). It 
increases with increasing coverage c at low c and decreases at higher c > 0.42. On the other 
hand, the mean energy associated with the vertical motion probably declines with increasing 
c for all c due to additional dissipation from collisions. In this interpretation the decrease 
in T with c at high c mirrors the decrease in energy associated with vertical motion. In 
contrast, the smaller values of T found at low c likely indicate that the energy transfer from 
vertical to horizontal motion becomes less efficient at low c < 0.42. 

We have shown that particles of different mass do not have the same kinetic energy or 
the same granular temperature when both are present simultaneously (Fig. |IB|), an apparent 
violation of equipartition. The most reasonable explanation is that all particles acquire 
similar vertical velocity fluctuations from the container. The more massive particles therefore 
obtain a larger vertical kinetic energy. This excess vertical energy is transferred to the 
horizontal motion, leading to a violation of equipartition. 

An important result of this investigation is that in the steady state the velocity dis- 
tributions deviate measurably from a Gaussian (Fig. |7|), but can be described well by 
P(v) ~ exp(— \vfvcl 1 - 5 ) for broad ranges of frequencies /, accelerations a and coverages 
c, in agreement with the theory of Ref ||. In most of our experiments the time between 
interparticle collisions r c is comparable to the time between contacts with the plate. How- 
ever, if forcing collisions become significantly less frequent than interparticle collisions (as for 



free cooling, Fig. [II]), then the distribution approaches an exponential. This experimental 
observation is consistent with the numerical results of Puglisi et al. [[F]. This quantitative 
correspondence with experiment indicates that the theoretical and numerical approaches of 
Refs. 0,0] to describing the statistical properties of granular particles are promising. 

Puglisi et al. suggest that non-Gaussian behavior and clustering are indications of 
essentially the same particle dynamics. If so, our experimental results indicate that clustering 
must occur for a very large range of excitation amplitudes and frequencies. 
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TABLES 



Material 


Diameter 


Mass 


Density 




[cm] 


[g] 


[g/cm 3 ] 


Glass (black and white) 


0.407 


0.0864 


2.46 


Stainless Steel 316 


0.159 


0.0168 


8.03 


Stainless Steel 316 


0.318 


0.1345 


8.03 


Stainless Steel 302 


0.476 


0.4539 


8.03 



TABLE I. Material parameters of granular particles. Black glass beads are used as background 
particles in all experiments. 
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FIGURES 



FIG. 1. Measured trajectories of some white particles moving among a sea of dark ones at a 
total coverage c = 0.28. The container is vertically vibrated at 100 Hz and peak acceleration 5 g. 
We acquire 546 frames at 250 frames/s. The particle diameter is about 11 pixels). 

FIG. 2. Apparent non-dimensional granular temperature T(At) corresponding to displacements 
along one axis in the sample interval as a function of sampling rate 1/At. The approach to 
saturation for large sampling rates occurs when the sampling rate far exceeds the collision rate. In 
the opposite limit T decreases strongly. The lines are explained in the text. 

FIG. 3. Apparent non-dimensional granular temperature T(At) corresponding to displacements 
along one axis in the sample interval At as a function of sampling rate 1 /At for a = 5 g. The 
limiting value for large sampling rates depends on coverage. 

FIG. 4. Granular temperature T determined from displacements along one axis as a function 
of coverage for a range of peak plate accelerations. 

FIG. 5. Non-dimensional granular temperature T vs. coverage. It is approximately indepen- 
dent of acceleration. 

FIG. 6. Granular temperature T vs. vibration frequency / at three coverages. It declines 
approximately as l// 2 . 

FIG. 7. Velocity distribution of the tracked particles plotted on linear (a,c) and logarithmic 
(b,d) scales (driving acceleration a = 5 g, coverage c = 0.42, and frequency / = 100 Hz). The 
solid line is a fit to Fi (Eq. ||). The Gaussian fit (dashed line) underestimates the probablilities at 
high velocities. The data were weighted equally on linear scales in (a,b) and equally on logarithmic 
scales in (c,d), respectively. 

FIG. 8. Distributions of non-dimensional velocity v obtained from displacements along one 
direction vs. v for (a) a range of accelerations and (b) a range of coverages. Fits to Eq. I (1 mes) 
are also shown. Data are shifted vertically in some cases for clarity. Deviations from Eq. |5| occur 
for c < 0.28. 

FIG. 9. Distributions of non-dimensional velocity v obtained from displacements along one 
direction vs. v 1 ' 5 at different vibration frequencies, and fits to Eq. ||. (c = 0.42, a = 5g) 

FIG. 10. Non-dimensional velocity distribution at low acceleration a = 2 g, and fits to 
P(v) ~ exp{— \v\ a ) showing the crossover from a Gaussian to an exponential function. 
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FIG. 11. Time averaged velocity distribution during free cooling after excitation at a = 5 g 
and / = 100 Hz (c = 0.84). The high velocity tail is exponential (a straight line on this log- linear 
plot). 

FIG. 12. Velocity distributions for steel beads moving among glass beads (c = 0.42, a = 4 g, 
/ = 100 Hz). Distributions are fitted to Eq. [| (solid lines) and also to a Gaussian (dashed line) for 
the largest beads. 

FIG. 13. Mean kinetic energy of steel beads moving among glass beads (scaled by the energy 
of the surrounding glass beads) as a function of particle mass, (a = 4 g; / = 100 Hz; squares, 
c = 0.42; triangles c = 0.84) 
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FIG. 1. Measured trajectories of some white particles moving among a sea of dark ones at a 
total coverage c = 0.28. The container is vertically vibrated at 100 Hz and peak acceleration 5 g. 
We acquire 546 frames at 250 frames/s. The particle diameter is about 11 pixels). 
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FIG. 2. Apparent non-dimensional granular temperature T(At) corresponding to displacements 
along one axis in the sample interval as a function of sampling rate 1/At. The approach to 
saturation for large sampling rates occurs when the sampling rate far exceeds the collision rate. In 
the opposite limit T decreases strongly. The lines are explained in the text. 
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FIG. 3. Apparent non-dimensional granular temperature T(At) corresponding to displacements 
along one axis in the sample interval At as a function of sampling rate I /At for a = 5 g. The 
limiting value for large sampling rates depends on coverage. 
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FIG. 5. Non-dimensional granular temperature T vs. coverage. It is approximately indepen- 
dent of acceleration. 
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FIG. 7. Velocity distribution of the tracked particles plotted on linear (a,c) and logarithmic 
(b,d) scales (driving acceleration a = 5 g, coverage c = 0.42, and frequency / = 100 Hz). The 
solid line is a fit to F2 (Eq. ||). The Gaussian fit (dashed line) underestimates the probablilities at 
high velocities. The data were weighted equally on linear scales in (a,b) and equally on logarithmic 
scales in (c,d), respectively. 
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FIG. 8. Distributions of non-dimensional velocity v obtained from displacements along one 
direction vs. v for (a) a range of accelerations and (b) a range of coverages. Fits to Eq. | (1 mes) 
are also shown. Data are shifted vertically in some cases for clarity. Deviations from Eq. |5| occur 
for c < 0.28. 
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FIG. 10. Non-dimensional velocity distribution at low acceleration a = 2 g, and fits to 
P(v) ~ exp{— \v\ a ) showing the crossover from a Gaussian to an exponential function. 
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FIG. 11. Time averaged velocity distribution during free cooling after excitation at a = 5 g 
and / = 100 Hz (c = 0.84). The high velocity tail is exponential (a straight line on this log- linear 
plot). 
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FIG. 12. Velocity distributions for steel beads moving among glass beads (c = 0.42, a = 4 g, 
/ = 100 Hz). Distributions are fitted to Eq. [| (solid lines) and also to a Gaussian (dashed line) for 
the largest beads. 



26 




0.1 



10 



m stee/ m glass 



FIG. 13. Mean kinetic energy of steel beads moving among glass beads (scaled by the energy 
of the surrounding glass beads) as a function of particle mass, (a = 4 g; / = 100 Hz; squares, 
c = 0.42; triangles c = 0.84) 
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